Spurious diapycnal mixing in terrain-following coordinate models: The problem and a solution   总被引：5，自引：0，他引：5  

Patrick Marchesiello  Laurent Debreu  Xavier Couvelard 《Ocean Modelling》2009,26(3-4):156-169

In this paper, we identify a crucial numerical problem in sigma coordinate models, leading to unacceptable spurious diapycnal mixing. This error is a by-product of recent advances in numerical methods, namely the implementation of high-order diffusive advection schemes. In the case of ROMS, spurious mixing is produced by its third-order upwind advection scheme, but our analysis suggests that all diffusive advection schemes would behave similarly in all sigma models. We show that the common idea that spurious mixing decreases with resolution is generally false. In a coarse-resolution regime, spurious mixing increases as resolution is refined, and may reach its peak value when eddy-driven lateral mixing becomes explicitly resolved. At finer resolution, diffusivities are expected to decrease but with values that only become acceptable at resolutions finer than the kilometer. The solution to this problem requires a specifically designed advection scheme. We propose and validate the RSUP3 scheme, where diffusion is split from advection and is represented by a rotated biharmonic diffusion scheme with flow-dependent hyperdiffusivity satisfying the Peclet constraint. The rotated diffusion operator is designed for numerical stability, which includes improvements of linear stability limits and a clipping method adapted to the sigma-coordinate. Realistic model experiments in a southwest Pacific configuration show that RSUP3 is able to preserve low dispersion and diffusion capabilities of the original third-order upwind scheme, while preserving water mass characteristics. There are residual errors from the rotated diffusion operator, but they remain acceptable. The use of a constant diffusivity rather than the Peclet hyperdiffusivity tends to increase these residual errors which become unacceptable with Laplacian diffusion. Finally, we have left some options open concerning the use of time filters as an alternative to spatial diffusion. A temporal discretization approach to the present problem (including implicit discretization) will be reported in a following paper.  相似文献   

Mechanics of material interfaces,edited by A. P. S. Selvadurai and G. Z. Voyiadjis,Studies in Applied Mechanics No. 11, Elsevier Amsterdam, 1986. No. of pages: 245. ISBN 0-444-42625-6 (Vol. 11)     

Laurent Vulliet 《国际地质力学数值与分析法杂志》1989,13(1):108-108

  相似文献   

Outgassing-induced effects in the rotational state of comet 67P/Churyumov—Gerasimenko during the Rosetta mission     

Pedro J. Gutiérrez  Laurent Jorda  Philippe Lamy 《Planetary and Space Science》2005,53(11):1135-1145

The new target of the Rosetta mission is comet 67P/Churyumov-Gerasimenko (hereafter 67P/C-G). In order to support the planning of the mission, in particular the strategy during the mapping and landing phases, we have performed numerical simulations of the rotational evolution of a comet in the orbit of 67P/C-G. In these simulations, the currently known observational constraints have been taken into account and a large set of initial conditions were considered. For most of the simulations, we observe that the sublimation-induced torques produce significant changes in the rotational parameters of a 67P/C-G-like comet. Typical rates of change for the spin period from the rendezvous up to the end of the nominal mission range from 0.001 to depending on different circumstances as described in the text. At perihelion, rates of change of the orientation of the angular momentum vector amount to about 0.01-. These simulations suggest that a specific strategy should be defined in order to monitor likely variations of the rotational parameters. As an example we show a possible optimized schedule for observations with the OSIRIS instrument to determine the rotational parameters of comet 67P/C-G and their possible evolution.  相似文献   

The effect of Congo River freshwater discharge on Eastern Equatorial Atlantic climate variability     

Stefano Materia  Silvio Gualdi  Antonio Navarra  Laurent Terray 《Climate Dynamics》2012,39(9-10):2109-2125

The surface ocean explains a considerable part of the inter-annual Tropical Atlantic variability. The present work makes use of observational datasets to investigate the effect of freshwater flow on sea surface salinity (SSS) and temperature (SST) in the Gulf of Guinea. In particular, the Congo River discharges a huge amount of freshwater into the ocean, affecting SSS in the Eastern Equatorial Atlantic (EEA) and stratifying the surface layers. The hypothesis is that an excess of river runoff emphasize stratification, influencing the ocean temperature. In fact, our findings show that SSTs in the Gulf of Guinea are warmer in summers following an anomalously high Congo spring discharge. Vice versa, when the river discharges low freshwater, a cold anomaly appears in the Gulf. The response of SST is not linear: temperature anomalies are considerable and long-lasting in the event of large freshwater flow, while in dry years they are less remarkable, although still significant. An excess of freshwater seems able to form a barrier layer, which inhibits vertical mixing and the entrainment of the cold thermocline water into the surface. Other processes may contribute to SST variability, among which the net input of atmospheric freshwater falling over EEA. Likewise the case of continental runoff from Congo River, warm anomalies occur after anomalously rainy seasons and low temperatures follow dry seasons, confirming the effect of freshwater on SST. However, the two sources of freshwater anomaly are not in phase, so that it is possible to split between atypical SST following continental freshwater anomalies and rainfall anomalies. Also, variations in air-sea fluxes can produce heating and cooling of the Gulf of Guinea. Nevertheless, atypical SSTs cannot be ascribed to fluxes, since the temperature variation induced by them is not sufficient to explain the SST anomalies appearing in the Gulf after anomalous peak discharges. The interaction processes between river runoff, sea surface salinity and temperature play an effective role in the interannual variability in the EEA region. Our results add a new source of variability in the area, which was often neglected by previous studies.  相似文献   

Strategies to mitigate aliasing of loading signals while estimating GPS frame parameters   总被引：7，自引：5，他引：2  

Xavier?CollilieuxEmail author  Tonie?van?Dam  Jim?Ray  David?Coulot  Laurent?Métivier  Zuheir?Altamimi 《Journal of Geodesy》2012,86(1):1-14

Although GNSS techniques are theoretically sensitive to the Earth center of mass, it is often preferable to remove intrinsic origin and scale information from the estimated station positions since they are known to be affected by systematic errors. This is usually done by estimating the parameters of a linearized similarity transformation which relates the quasi-instantaneous frames to a long-term frame such as the International Terrestrial Reference Frame (ITRF). It is well known that non-linear station motions can partially alias into these parameters. We discuss in this paper some procedures that may allow reducing these aliasing effects in the case of the GPS techniques. The options include the use of well-distributed sub-networks for the frame transformation estimation, the use of site loading corrections, a modification of the stochastic model by downweighting heights, or the joint estimation of the low degrees of the deformation field. We confirm that the standard approach consisting of estimating the transformation over the whole network is particularly harmful for the loading signals if the network is not well distributed. Downweighting the height component, using a uniform sub-network, or estimating the deformation field perform similarly in drastically reducing the amplitude of the aliasing effect. The application of these methods to reprocessed GPS terrestrial frames permits an assessment of the level of agreement between GPS and our loading model, which is found to be about 1.5 mm WRMS in height and 0.8 mm WRMS in the horizontal at the annual frequency. Aliased loading signals are not the main source of discrepancies between loading displacement models and GPS position time series.  相似文献   

Uncertainties in European summer precipitation changes: role of large scale circulation   总被引：2，自引：2，他引：0  

Julien Boé  Laurent Terray  Christophe Cassou  Julien Najac 《Climate Dynamics》2009,33(2-3):265-276

Climate models suggest that anthropogenic emissions are likely to induce an important drying during summer over most of Europe in the late 21st century. However, the amplitude of the associated decrease in precipitation strongly varies among the different climate models. In order to reduce this spread, it is first necessary to identify its causes and the associated physical mechanisms. Consequently, the focus of this paper is to better estimate the role of large scale circulation (LSC) in precipitation changes over Europe using a multi-model framework and then to characterize the LSC changes using the weather regime paradigm. We show that LSC changes directly lead to a decrease of precipitation over northwestern Europe. This circulation-driven decrease in rainfall is mainly linked to an increase (decrease) of the occurrence of positive (negative) phase of the North Atlantic Oscillation regime. LSC is also responsible for a significant part of the models spread in precipitation changes over these regions. Over southern Europe, the role of LSC changes on multi-model mean precipitation changes is generally weak. We also show that the precipitation anomalies directly induced by LSC modifications seem to be further amplified through local feedbacks.  相似文献   

A chemical and thermodynamic model of oil generation in hydrocarbon source rocks     

Harold C. Helgeson  Laurent Richard  William F. McKenzie  Alexandra Schmitt 《Geochimica et cosmochimica acta》2009,73(3):594-7350

Thermodynamic calculations and Gibbs free energy minimization computer experiments strongly support the hypothesis that kerogen maturation and oil generation are inevitable consequences of oxidation/reduction disproportionation reactions caused by prograde metamorphism of hydrocarbon source rocks with increasing depth of burial.These experiments indicate that oxygen and hydrogen are conserved in the process.Accordingly, if water is stable and present in the source rock at temperatures ?25 but ?100 °C along a typical US Gulf Coast geotherm, immature (reduced) kerogen with a given atomic hydrogen to carbon ratio (H/C) melts incongruently with increasing temperature and depth of burial to produce a metastable equilibrium phase assemblage consisting of naphthenic/biomarker-rich crude oil, a type-II/III kerogen with an atomic hydrogen/carbon ratio (H/C) of ∼1, and water. Hence, this incongruent melting process promotes diagenetic reaction of detritus in the source rock to form authigenic mineral assemblages.However, in the water-absent region of the system CHO (which is extensive), any water initially present or subsequently entering the source rock is consumed by reaction with the most mature kerogen with the lowest H/C it encounters to form CO₂ gas and a new kerogen with higher H/C and O/C, both of which are in metastable equilibrium with one another.This hydrolytic disproportionation process progressively increases both the concentration of the solute in the aqueous phase, and the oil generation potential of the source rock; i.e., the new kerogen can then produce more crude oil.Petroleum is generated with increasing temperature and depth of burial of hydrocarbon source rocks in which water is not stable in the system CHO by a series of irreversible disproportionation reactions in which kerogens with higher (H/C)s melt incongruently to produce metastable equilibrium assemblages consisting of crude oil, CO₂ gas, and a more mature (oxidized) kerogen with a lower H/C which in turn melts incongruently with further burial to produce more crude oil, CO₂ gas, and a kerogen with a lower H/C and so forth.The petroleum generated in the process progresses from heavy naphthenic crude oils at low temperatures to mature petroleum at ∼150 °C. For example, the results of Computer Experiment 27 (see below) indicate that the overall incongruent melting reaction in the water-absent region of the system C-H-O at 150 °C and a depth of ∼4.3 km of an immature type-II/III kerogen with a bulk composition represented by C₂₉₂H₂₈₈O_12(c) to produce a mature (oxidized) kerogen represented by C₁₂₈H₆₈O_7(c), together with a typical crude oil with an average metastable equilibrium composition corresponding to C_8.8H_16.9 (C_8.8H_16.9(l)) and CO₂ gas (CO_2(g)) can be described by writing